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Abstract 

The Relativistic Mean Field (RMF) approach describing the motion of independent particles in 
effective meson fields is extended by a microscopic theory of particle vibrational coupling. It leads 
to an energy dependence of the relativistic mass operator in the Dyson equation for the single- 
particle propagator. This equation is solved in the shell-model of Dirac states. As a result of the 
dynamics of particle- vibrational coupling we observe a noticeable increase of the level density near 
the Fermi surface. The shifts of the single-particle levels in the odd nuclei surrounding 208 Pb and 
the corresponding distributions of the single-particle strength are discussed and compared with 
experimental data. 

PACS numbers: 21.30.Ee, 21.60.Jz, 21.65. +f, 21.10.-k 



I. INTRODUCTION 



New experimental facilities with radioactive nuclear beams make it possible to investigate 
the nuclear chart not only along the narrow line of stable isotopes but also in areas of large 
neutron- and proton excess far from the valley /^-stability. This situation has stimulated 
enhanced efforts on the theoretical side to understand the dynamics of the nuclear many- 
body problem by microscopic methods. Very light nuclei with A < 12 are studied by an "ab 
initio" approach utilizing bare nucleon-nucleon interactions of two- and three-body character 
and modern shell-model calculations based on large scale diagonalization techniques and 
truncation schemes show considerable success in predicting properties of somewhat heavier 
nuclei. For the large majority of nuclei, however, a quantitative microscopic description is 
only possible by density functional theory. These methods are based on mean-field theory. 
Although density functional theory can, in principle, provide an exact description of the 
many-body dynamics, if the exact density functional is known, in nuclear physics one is 
far from a microscopic derivation of this functional. The most successful schemes use a 
phenomenological ansatz incorporating as many symmetries of the system as possible and 
adjust the parameters of these functionals to ground state properties of characteristic nuclei 
all over the periodic table. Considerable progress has been reported recently in constructing 
such functionals. For a recent review see [jj. 

Of particular interest are covariant density functionals because they are based on 

Lorentz invariance. This symmetry not only allows to describe the spin-orbit coupling, 
which has an essential influence on the underlying shell structure, in a consistent way, but it 
also put stringent considerable restrictions on the number of parameters in the corresponding 
functionals without reducing the quality of the agreement with experimental data. A very 
successful example is the Relativistic Hartree-Bogoliubov model |4|. It combines a density 
dependence through a non-linear coupling between the meson fields [5] with pairing correla- 
tions based on an effective interaction of finite range. A large variety of nuclear phenomena 
have been described over the years within this model: the equation of state in symmetric 
nuclear matter,ground state properties of finite spherical and deformed nuclei all over the 
periodic table ff from light nuclei Q to super-heavy elements g from the neutron dnp 
line, where halo phenomena are observed [2[ to the proton drip line |H| with nuclei unsta- 
ble against the emission of protons [ljj. Recently this model has been also applied very 
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successfully for the de scrip tion of excited states, such as rotational bands in normal and 
super-deformed nuclei |l2i ll^ and collective vibrations jl^. Rotations are treated in the 
cranking approximation, which provides a quasi-static description of the nuclear dynamics 
in a rotating frame and for the description of vibrations a time-dependent mean field approx- 
imation is used by assuming independent particle motion in time-dependent average fields 
[15^ . In the small amplitude limit one obtains the relativistic Random Phase Approximation 
(RRPA) jl6l| . This method provides a natural framework to investigate collective and non- 
collective excitations of p/i-character. It is successful in particular for the understanding of 
the position of giant resonances and spin- or/and isospin-excitations as the Gamov Teller 
Resonance (GTR) or the Isobaric Analog Resonance (IAR). Recently it has been also used 
for a theoretical interpretation of low-lying El-strengths observed in neutron rich isotopes 
(pygmy modes) [l^l and for low-lying collective quadrupole excitations 
Of course the density functional theory based on the mean field framework cannot provide 
an exact treatment of the full nuclear dynamics. It is known to break down not only in 
transitional nuclei, where one has to include correlations going beyond the mean field ap- 
proximation by treating quantum fluctuations through a superposition of several mean field 
solutions, as for instance in the Generator Coordinate Method (GCM) 19], but already in 
ideal shell model nuclei such as 208 Pb with closed protons and neutron shells one finds in 
self-consistent mean field calculations usually a single particle spectrum with a considerably 
enhanced Hartree-Fock gap in the spectrum and a reduced level density at the Fermi surface 
as compared with the experiment. It is well known that this fact is connected with the 
relatively small effective mass in such models. 

n 

Mahaux and collaborators [20] have shown that the effective mass in nuclear matter is 
roughly m*/m pa 0.8. In finite nuclei it should be modified by the coupling of the single 
particle motion to low-lying collective surface vibrations. This leads, in the vicinity of the 
Fermi surface, to an enhancement of m*/m pa 1. Non-self-consistent models with the bare 
mass ( m*/m pa 1) show indeed a single particle spectrum with a level density close to the 
experiment. 

Using the quasi-particle concept of Landau theory and Green's function techniques, one can 
derive a one-body equation for the single-particle Green's function, which is in principle 



exact, the Dyson equation 



2l\ . It contains a non-local and energy dependent self-energy, 



also called mass-operator. The energy independent part of this operator can be described 
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very well in mean field theory. The most important origin of an energy dependence is given 
by the coupling of the single particle motion to low-lying collective vibrations. 



II. THE ENERGY DEPENDENT PART OF THE MASS OPERATOR 
A. One-nucleon motion in the relativistic mean field 

In the relativistic nuclear theory the motion of the nucleons is described by the Dirac equation 

(7^-7^)1^ = 0, (1) 

where the effective mass is given by 

m* = m + S s (2) 

with the scalar part S s of the mass operator and where the generalized four- vector momentum 
operator has the form 

p^ = p ^-j:^= (i__£ 0j iv + £j (3) 

with the vector part S M of the mass operator 

£" = (£ ,£). (4) 

The index V in the Eq. (J2J) denotes that the effective mass is described by the scalar o- 
meson field. In order to characterize ground state properties the stationary Dirac equation 
has to be solved: 

(a(p-£)+/3m* + £ )|V>>=<#>. (5) 
In the general case the full mass operator is non-local in the space coordinates and also in 
time. This non-locality means that its Fourier transform has both momentum and energy 
dependence. Let us assume the components of mean field to be sums of the stationary local 
and energy dependent non-local terms: 

S(r, r'; u) = S(r)5(r - r') + S e (r, r'; u), (6) 

where all the components of the mass operator are involved: 

£ = (£„£") 
£ = (£ S ,£ M ) 

£ e = (£^,£^) (7) 
4 



and index "e" indicates the energy dependence. 

The scalar component of the energy-independent mass operator is proportional to the o~- 
meson field: 

t s (r)=g a a(r). (8) 

Time-like and space-like components of the local and energy-independent part of the mass 
operator are generated by the isoscalar w-meson and isovector p-meson fields co> M ,p M 
and Coulomb field 

W(r) =g^(r)+g p Tp^r) + e^-^-A^T), (9) 

where arrows denote isovectors and bold-faced letters indicate vectors in three-dimensional 
space. These fields satisfy the inhomogeneous Klein-Gordon equations: 

(-A + ml) a(r) = -g aPs (r) - g 2 a 2 (r) - g 3 a 3 (r) (10) 

(-A + m^)^(r) = ^/(r) (11) 

(-A + m 2 p ) p»(r) = g p j»(r) (12) 

-A#(r) = e#(r), (13) 

where the sources are determined by the respective density and current distributions in a 
system of A nucleons: the scalar density for a-field 

A 

p s (r) = ^^W^(r), (14) 
i=i 

the baryon current for the cu-field 

A 

^(r)=^^(r) 7 ^(r), (15) 



i=i 



the isovector current for the p-field 
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;r) = ]T^(r)7^(r), (16) 
i=i 

and the charge current for the photon-field 

j;W = E^(r)7 /i ^^^(r). (17) 
i=i 

The summation in (|14I-I16|) is performed over occupied states in the Fermi sea, in accordance 
with no-sea approximation, so that the contribution of the negative-energy states to the 
densities and currents is neglected. 



B. The single-particle Green's function 



In the present work we assume time-reversal symmetry that means the absence of currents 
in the nucleus and, thus, vanishing space-like components of E. The equation of the one- 
nucleon motion has the form: 

(h v + PZ(e) + 1%(e))\tl>)=e\1>) (18) 

or, in the language of Green's functions 

(e-h v -WZ(e)-I%(e))G(e) = l, (19) 

where h v denotes the Dirac hamiltonian with the energy-independent mean field: 

h v = ap + (3{m + E s ) + E . (20) 

We now work in the shell-model Dirac basis {l^fc)} which diagonalizes the energy- 
independent part of the Dirac equation: 

h v \rl> k )=e k \i(> k ). (21) 

In addition we assume in the present work spherical symmetry. In this case the spinor 
\ipk) is characterized by the set of single-particle quantum numbers k = {(k),m k },(k) = 
{t k , ilk, rik,jk, k} with the radial quantum number n k , angular momentum quantum numbers 
jk,m k , parity n k and isospin t k : 

if ,\ ( fm(r)Qi kjkmk (&,<p,s) \ 

M T >s,t) = , , , « J Xt k (t), (22) 

where the orbital angular momenta l k and l k of the large and small components are deter- 
mined by the parity of the state k: 

k = 3k + \, k =3k~\ for n k = (-l) jk+ ^ 
k = jk - \, k = jk + \ for ir k = (-l) Jfc_ 5, 

/(fe)(r) and g( k )(r) are radial wave functions and <&y m is a two-dimensional spinor: 

$ ljm (0,(p,s)= ^2{ 1 -m s lmi\jm)Y lm[ (^,(f)Xm s (s)- (24) 

m s mi 
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In this basis one can rewrite Eq. (|19j) as follows: 

J2H £ - - £|,(e)}GW(e) = hk>, (25) 
i 

where the letter indices k, k', I denote full sets of the spherical quantum numbers mentioned 
above. 

In the next step we represent the exact single-particle Green's function entering Eq. t)19j) 
in the Lehmann expansion. In contrast to the non-relativistic case, where one has occupied 
states below the Fermi surface (hole states h) and empty states above the Fermi surface 
(particle states p) we now have according to the no-sea approximation in addition empty 
states with negative energies (antiparticle states a). For a detailed discussion of this point 



see also Ref. 
form 



|. Therefore the Lehmann representation of the Green's function has the 

,A0, AO* Op Op* ,.0a,.0a* 

Gkl ( £ ) = y J^li — + y Jihli — + y J^l — (26 ) 

L - J e - e h - in / — ' e - e„ + in ^— ' e - e a + in 

h p r a 

where n — > +0 and the matrix elements are defined as 

Xt° = (h\^ k \0), 



xl p = m k \ P ), 

x l a = (0\^ k \a). (27) 

Here ipk is the Dirac field annihilation operator of the state k. |0) denotes the ground state 
of the subsystem of N particles in the even-even nucleus in no-sea approximation, i.e. the 
negative energy states are essentially empty. The states \h) correspond to the ground state 
and to excited states of the subsystem of (TV — 1) particles and \p) are the ground and excited 
states of the system of (N + 1) particles, respectively. Because of the no-sea approximation 
the negative energy states \a) are not occupied in |0) and therefore there exist also states 
| a) in the (N + 1) particle system where a level with negative energy are occupied. 
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C. The pole structure of the mass operator 



Let us now define the energy-dependent part of the mass operator (simply called 'mass 
operator' in the following). Its matrix elements have the form: 

ES,(e) = J M'^(r)(^(r,r'; £ )+^(r,r'i£))^(r'), (28) 

Obviously, on this stage one needs some model assumptions. In the present work we choose 
a rather simple particle-phonon coupling model 22] to describe the energy dependence of 
S e . Within this model S e is a convolution of the particle-phonon coupling amplitude T and 
the exact single-particle Green's function j^j: 



-oo 



/duo 
—r k iti k i(u})G k 'i>(e + uj), (29) 
Z711 




k'l 

— oc 



where the amplitude T has the following spectral expansion 

rmk , (u) = -Y( -J^l 7kkLtSL_) (30) 

n 

in terms of phonon vertexes 7^ and their frequencies Q^. They are determined by the 
following relation: 

Ha = ^2 v MwSp^ v . (31) 

k'l' 

Vki'ik' denotes the relativistic matrix element of the residual interaction and Sp is the transi- 
tion density. In the present work we use the linearized version of the model which assumes 
that 5p is not influenced by the particle-phonon coupling and can be computed within rela- 
tivistic RPA. The linearized version implies also that the energy- dependent part of the mass 
operator (J28|) contains the mean field Green's function G(e) = (e — h®)^ 1 instead of the 
exact Green's function G. So, the equation (|25|) becomes linear with respect to G. Since the 
mean field Green's function is 

G k i{e) = , (32) 

where a k = +1 if k is an unoccupied state of p- or a-types and a k = — 1 for an occupied k 
states of /i-type, the mass operator S e takes the form: 

ye ( s *,•,;>,/■ (k || 7^" } II n)(l ]| 7^_j II n)* 
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FIG. 1: The particle £y „ and the hole S^ /?i „ components of the relativistic mass operator in the 
graphical representation, p, a, h are the particle, antiparticle and hole types of the intermediate 
states. See text for the detailed explanation. 

Here we use the notation: 

(k || Y {an) II n) = 5 antl (k || Y II n) + d^in \\ Y II &)*■ (34) 

Since the indexes k, I and n run through the whole Dirac basis, each state in ()3*3*)) can be 
a particle above the Fermi surface, a hole below the Fermi surface or a particle in a state 
with negative energy (antiparticle state). The graphical representation of the mass operator 
is given in Fig. ^ We draw the particle and the hole components assuming all the possible 
types of intermediate states. Solid line with arrow denotes a particle (hole) in the Fermi sea, 
dashed line means a particle in the empty Dirac sea, weavy line is a phonon propagator, and 
small circle denotes a phonon vertex 1)31)1 . Time direction is from the left to the right. 
One can see from the Eq. that the matrix E|j contains a small number of the off- 

diagonal elements with relatively large energy denominators. Additionally, it was shown by 
explicit calculations within the non- relativistic approach |2J] that it is justified to use the 
diagonal approximation: 

E e w (e) = S^Ke) (35) 
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with 



1 ^ I (A; || 
. + 1 ^ e - e n - a n ({ 



n)| 5 



(36) 



2 Jk + 1 8 ~ £ n ~ ^n(^ ~ iff) 

In analogy with the conventional terminology of non-relativistic approaches, let us call the 
intermediate term n 'polarization term' if a n = and 'correlation term' if a n = —<J k - The 
correlation term describes, obviously, the backwards going diagrams in Feynman's language 
and corresponds to the ground state correlations caused by the particle-vibration coupling. 
Thus, within the diagonal approximation of the mass operator ()35j) the exact Green's function 
G is also diagonal in the Dirac basis and the Dyson equation forms for each k a non-linear 
eigenvalue equation 

(e-e fc -E|(e))G fc (e) = l. (37) 
The poles of the Green's function Gk{s) correspond to the zeros of the function 



f(e) = e-e k -Zl 



(3* 



For each quantum number k there exist several solutions characterized by the index A. 
Because of the coupling to the collective vibrations the single particle state k is fragmented. 
In the vicinity of the pole the Green's function can be represented as follows: 



G 



(A), 



s, 



(A) 



e — e 



(A) 



+ icTkV 



(39) 



where the residuum has a meaning of the single-particle (hole) strength of the state A 
with single-particle quantum numbers k. Differentiation of the equation (|37j) with respect 
to e provides the expression for the residua: 

-l 



S 



(A) 



1 



de 



(40) 



There are several ways to solve the equation (J37|) . In the present work we employ the 



method which has been used in Ref. 



24j to solve the similar problem in the non-relativistic 



framework. Since the mass operator of the form ()33|) has a simple-pole structure, it is 
convenient to reduce the Eq. (}3*Tj) to a diagonalization problem of the following matrix: 



£ k 



'lkn 2 



\ 



CI + e ni 

a n2 n^+e n2 




Ml * 
Vkn 2 








(41) 
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where 

(k II || m) , , 

The eigenvalues of the matrix (|41|) are the desired poles ej^ of the exact Green's function. 
The structure of the solution is well known: these eigenvalues lie between the poles of the 
mass operator. Eventually, the spectroscopic factors have to be calculated at the points of 
these poles according to (J4"U|) : 

sw = (i+ 1 ^ 1(^11 7^ INI 2 



— V Iffll^'lln)!' V\ (43) 



2 ^ + 1 ^ (4 J - ^ - 

III. DETAILS OF THE CALCULATIONS AND DISCUSSION 

The matrices (j4Tj) have been diagonalized for the single-particle states ki of both neutron 
and proton subsystems belonging to the four major shells around iV = 126 and Z = 82. 
Thus, the eigenvalues with the largest spectroscopic factors correspond to the single-particle 
excitations of the nuclei 207 Pb, 209 Pb, 207 T1 and 209 Bi. In subsection 1111 Al we discuss the 
effect of states with negative energies in the Dirac sea on the mass operator relying on results 
obtained within the restricted particle-phonon space. More realistic results for energies 
and spectroscopic factors obtained in an enlarged particle-phonon space are presented and 
discussed in subsection IIII Bl In subsection IHI CI we compare the present method with other 
approaches. 



A. Relativistic effects: illustrative calculations 



The main interest of the present work is to describe the effects of complex configurations 
within the relativistic scheme. Therefore, first we investigate the contributions of pure 
relativistic terms to the mass operator and, hence, their influence on the single-particle 
spectrum of odd nuclei. In order to keep the numerical effort within a reasonable limit we 
used a restricted particle-phonon space taking into account only the most collective phonons 
with spin and parity J 71 = 2 + , 3~, 4 + , 5~, 6 + below the neutron separation energy and 
a reduced number of single particle states with positive energy (particles or holes). This 
enables one to reduce strongly the number of poles in the mass operator (fHHJ) as well as the 
dimension of the matrix (|41jl . Notice, that since in the Green's function formalism we stay in 
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the single-particle basis, it is always possible to vary the number of phonons, and the problem 
of the completeness of the phonon basis does not arise at all. In all these calculations we 



25| for the Lagrangian. 



use the parameter set NL3 
The numerical results obtained in these investigations are compiled in the Table |l] For 
the first shell of neutron levels above ('particle') and below ('hole') the Fermi level three 
versions are given: in the version pha the index n in Eq. (J36j) includes all contributions 
from intermediate states above the Fermi level (p, with a n = +1), below the Fermi level (h, 
with cr n = —1) and in the Dirac sea (a, with a n = +1). Version pa (for particles) or h (for 
holes) excludes the backward going diagrams (i.e. only states with a n = are taken into 
account), and the third version ph does not contain antiparticle intermediate states in ()36|) . 
In this way, one can see that the effects of ground state correlations (GSC) caused by the 
particle-phonon coupling and neglected in the second version are significant and it is essential 
to take them into account in a realistic calculation. On the other hand, the contribution 
of the antiparticle subspace to the mass operator is quantitatively not of great importance. 
This can be understood by the fact that these configurations provide large values for the 
energy denominators in (|3*Bj) . Thus it is justified to disregard them in the full calculation. 
Notice, however, that version ph does not eliminate the effects of the Dirac sea completely 
since the phonon vertices still contain this contribution. As it has been discussed in Ref. 
Q these terras p,a y an aaportaat roie ia a proper —at of relativist* RPA. Otherwise 
it is not possible to obtain reasonable properties for the isoscalar modes within RRPA. 



B. Realistic calculations in an enlarged space 

In this section we neglect the effects of the Dirac sea, i.e. the intermediate index n in Eq. 
(J36J) runs only over particle states above the Fermi level and holes below the Fermi level. 
It has been found in section IHI Al that this is a very reasonable approximation. Since we 
do not have to include these contributions, we are now able to enlarge the particle-hole 
basis considerably by taking into account particle-hole configurations far away from the 
Fermi surface. In this case we increase the collectivity of the phonons and, consequently, 
the strength of the particle-vibrational coupling. The phonon basis was also enriched by 
including higher-lying modes up to 35 MeV, although these modes are not so important as 
the low-lying ones. Vibrations with the quantum numbers of spin and parity J n = 2 + , 3~, 
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TABLE I: Energies e k ' and spectroscopic factors S k ' of the dominant neutron levels in 208 Pb 
calculated in the strongly restricted particle-phonon space, pha denotes full the calculation, pa (h) 
is the version without backwards going terms, and ph is the version without contribution of the 
antiparticle states in the mass operator (see text for details). 



State k 


£k, MeV 




ef\ MeV 










Particle 




pha 


pa 


ph 


pha 


pa 


ph 


2g9/2 


-2.50 


-2.85 


-3.14 


-2.88 


0.89 


0.92 


0.89 


lill/2 


-2.97 


-2.82 


-3.20 


-2.90 


0.94 


0.97 


0.94 


1J15/2 


-0.48 


-1.16 


-1.33 


-1.21 


0.70 


0.74 


0.70 


3d5/2 


-0.63 


-0.96 


-1.05 


-0.98 


0.93 


0.94 


0.93 


4sl/2 


-0.36 


-0.88 


-0.92 


-0.89 


0.93 


0.93 


0.93 


2g7/2 


-0.56 


-0.71 


-0.90 


-0.76 


0.92 


0.94 


0.92 


3d3/2 


-0.02 


-0.35 


-0.42 


-0.37 


0.93 


0.93 


0.93 


Hole 




pha 


h 


ph 


pha 


h 


ph 


3pl/2 


-7.66 


-7.67 


-7.40 


-7.70 


0.96 


0.98 


0.96 


2f5/2 


-9.09 


-8.97 


-8.71 


-9.02 


0.93 


0.96 


0.93 


3p3/2 


-8.40 


-8.20 


-7.87 


-8.22 


0.90 


0.94 


0.90 


1H3/2 


-9.59 


-9.30 


-9.07 


-9.36 


0.90 


0.92 


0.89 


2f7/2 


-11.11 


-10.20 


-9.98 


-10.22 


0.72 


0.76 


0.72 


(lh9/2)i 


-13.38 


-13.32 


-13.23 


-13.34 


0.52 


0.47 


0.53 


(lh9/2) 2 




-12.48 


-12.42 


-12.49 


0.31 


0.39 


0.29 



4 + , 5~, 6 + were included in the phonon space. One should keep in mind, however, that in the 
solution of the RRPA-equations for the vibrational states besides the usual p/i-components 
a large number of a/i-components of the Dirac sea was included. Of course, as it is usually 
done in the RRPA calculations, both Fermi and Dirac subspaces were truncated at energies 
far away from the Fermi surface: in the present work we fix the limits e ph < 100 MeV and 
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£ah > —1800 MeV with respect to the positive continuum. The energies and B(EL)f values 
of the most collective phonon modes calculated with the parameter set NL3 are displayed in 
Table ITU together with some experimental data. 

As one can see from the Table |Hj the characteristics of low-lying modes obtained in RRPA 
with the parameter set NL3 are, in general, in accordance with experimental data. For 
the lowest 2 + ,4 + vibrations the B(EL)f values are in a good agreement with experimental 
ones, only the energies are slighty too high, whereas for the lowest 3~,5~ their energies are 
reproduced rather well and the B(E3)f value is to some extent overestimated, and the first 
6 + state is more collective within RRPA then the observed one. 

In the present calculations the phonon space was confined also by a criterion for the B(EL)f 
values: all modes with B(EL)f values less than 10% of the maximal one were neglected for 
2 + ,3~,4 + , and less than 20% for higher multipolarities. Nevertheless, the mass operator 
(J36j) has been calculated in a rather wide particle-phonon space and therefore the single- 
particle strength is distributed over many states. The typical dimension of the matrix (j41j) 
is about two thousand and it varies depending on the state k. As it was mentioned above, 
contributions of antiparticle states (n = a in the intermediate sum over n in the mass 
operator (}36|) ) were excluded because they provide large values of the energy denominators 
in (JHBJ)- On the other hand, the contributions of the correlation terms (i. e. the terms with 
a n = —<7k in the intermediate sum over n in the mass operator (J3H|) ) have been fully taken 
into account since they are found to be quantitatively important and they compensate to 
some extent the polarization terms. 

The final results of these calculations are compiled in the Table IIHI All the energies are 
related to the experimental ground states of the respective odd nuclei surrounding the doubly 
magic nucleus 208 Pb. The numbers denote the RMF single particle energies, s^ th are the 
eigenvalues of the matrix (J41j) with the maximal spectroscopic factor, and e^ exp are the 
experimentally observed excitation energies. We display here only the results for one major 
shell below and one shell above the Fermi surface because in the next shells almost all the 
single-particle levels turn out to be strongly fragmented due to phonon coupling and it is no 
longer possible to determine the dominant levels in these shells, in other words, the concept 
of Landau quasi-particles breaks down at energies far away from the Fermi level. 
The difference between Ek and e^ th is the shift of the single-particle level k caused by 
the coupling to collective surface vibrations. Notice, that almost all the levels are moving 
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TABLE II: Energies and reduced transition probabilities of the most collective vibrations in 208 Pb 
obtained within RRPA and experimental data from 





RRPA 






Exp. 


J 


uJ 


D(rjLi J | 


uJ 






^ivie v j 


\e im j 


^ivie v j 


^e im j 


2+ 


4.98 


2.69xl0 3 


4.07 


2.97xl0 3 




£ 0/1 

O.o4 


O.oZ X 1U 








O DO 

0.38 


1 oOwi r\3 








in /in 
1Z.4U 


4.UoX 1U 








zz.9b 


1 nO\ / 1 r\3 

l.Uox ID 






o — 



2.(4 


/ .4bX 1U 


cx 1 

z.bl 


5.40(30) X 11) 




4.yo 


O.ol X 1U 








7.29 


5.90xl0 4 








22.27 


6.12xl0 4 






A + 


4.yb 


1 on w 1 r\7 
l.oy x iu 


4.0Z 


1 on w 1 n7 
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TABLE III: Energies e$ and spectroscopic factors of the dominant single-particle levels in odd 



nuclei of the ^ u °Pb surroundings calculated in the wide particle-phonon space. The experimental 



data are taken from 
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downwards providing thus a considerably better agreement with experimental energies then 
the pure RMF states. One can see also from this table that the dominant neutron and proton 
levels obtained in these calculations have large spectroscopic factors and are therefore rather 
good single-particle states. This is in agreement with experiment. Nonetheless the single- 
particle strength is distributed over many levels. One can easily understand the origin of 
these large spectroscopic factors from the structure of solutions of the (|37|). Since each root 
lies between the two neighbouring poles of the mass operator, only one root can be found in 
the rather wide window 

£ (max) _ jy,(min) < ^ < £ {mm) + Qt i(mm) ^ ^ 

i.e -10.11 MeV < u < -1.20 MeV for neutron and -10.75 MeV < to < -1.06 MeV for proton 
subsystems, respectively. Therefore, for the single-particle state k near the Fermi surface 
other roots turn out to be far away because the lowest 3^ phonon has rather high energy 
in magic nuclei, therefore the mixing is weak and the respective spectroscopic factor (j43J) is 
close to one. On the contrary, if the state is near the limits or outside of the window, there 
are many other roots in the vicinity, and a strong mixing leads to the appearance of several 
levels with comparable strength. It can be easily understood why in open-shell nuclei such a 
mixing is much stronger near the Fermi surface: the window ()44|) is noticeably smaller due 
to both the smaller gap between the last occupied and the first unoccupied levels, and much 
the lower energy of the first 2 + phonons (see, for instance, 

Nevertheless, some dominant levels near the Fermi surface have noticeably reduced strength 
because of the particle-phonon coupling. This is also confirmed by experiment. One can 
see such a situation, for instance, for the lj'15/2 and the 2/7/2 neutron states in 209 Pb and 
207 Pb, respectively. Only for the state lh 9 / 2 in the 207 Pb we do not find good agreement: 
the spectroscopic factor is less then a half of the experimental value. One also should keep 
in mind that the experimental spectroscopic factors depend considerably on the parameters 
used in the DWBA analysis. The proton states are found to be somewhat more fragmented 
then the neutron states while the dominant levels for the protons above the Fermi surface 
are more strongly shifted relatively to the RMF values then in the those of the neutrons. 
The distributions of the single-particle strength for selected levels in 209 Pb, 209 Bi, 207 Pb, and 
207 T1 are represented in the Figs. El El IH El respectively. 

To illustrate the results of these calculations we chose one state of pronounced single-particle 
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FIG. 2: Single-particle strength distribution for the 3d 3 /2 (left panel) and IJ15/2 (right panel) states 
in 209 Pb obtained in the calculations (solid lines) and the experimental strengths of the respective 
dominant levels (dashed lines). 

nature (left panels) and one noticeably fragmented state (right panels) for each nucleus, both 
from the first major shell above and below the Fermi surface. As in Table UTTl all the energies 
are related to the ground state energy of the corresponding odd nucleus. As it was already 
mentioned, in the present calculations the single-particle strength is distributed over about 
two thousand states but most of them are vanishingly small, so only the states with the 
strength exceeding 10~ 3 are drown. The experimental strength of the dominant levels are 
shown with dashed lines. Some examples of the strongly fragmented states from the second 
major shells above and below the Fermi surface are shown in the Fig. |3 
To illustrate the shifts in the level schemes we show in Figs. E]and|B]the single-particle spectra 
for neutrons and protons. The spectra calculated with the energy-dependent correction 
(RMF+PVC) demonstrate a pronounced increase of the level density around the Fermi 
surface of 208 Pb both for neutron and proton subsystems comparatively the pure RMF 
spectra. In some cases it turned out to be possible to invert the order of levels and reproduce 
the observed sequence as one can see for the lj'15/2 and the 3d 5 / 2 neutron states. Another 
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FIG. 3: The same as in Fig. |2]but for the 2/ 7 / 2 (left panel) and H13/2 (right panel) states in 209 Bi. 
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FIG. 4: The same as in Fig. [2] but for the 1*13/2 (left panel) and 2/ 7 / 2 (right panel) states in 207 Pb. 
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FIG. 5: The same as in Fig. |2]but for the l/ixi/2 (left panel) and 2d 5 / 2 (right panel) states in 207 T1. 

and more important example is the inversion of the 2gg/2 and H11/2 neutron states (in the 
Fig. [7| they look coincided) which enables one to reproduce the spin of the 209 Pb ground 
state. 

In order to quantify these results we calculated the average distance between two levels in the 
spectrum shown in Figs. E|and|Hl (lh9/2 and lg7/2 states with small spectroscopic factors 
were excluded from the estimation of the neutron and proton spectrum, respectively). We 
obtain for the neutrons 1.0 (RMF), 0.83 (RMF+PCV) and 0.76 (EXP) in units of MeV. 
This corresponds to a level density of 1.0 (RMF), 1.20 (RMF+PCV) and 1.31 (EXP) in 
units of MeV -1 . The level density in the neighborhood of the Fermi surface is therefore in 
RMF-calculations by a factor 0.76 smaller than the experimental value. Taking into account 
particle-vibrational coupling we find only a reduction of 0.92. Assuming an effective mass 
close to 1 for the experiment, and taking into account that the level density at the Fermi 
surface is proportional to m*/m, this correspond to an effective mass m*/m ^ 0.76 for the 
RMF and m*/m « 0.92 for the RMF+PCV calculations. For the protons the situation is 
similar. From Fig. [Hjwe obtain for the average level distance 1.50 (RMF), 1.24 (RMF+PCV) 
and 1.06 (EXP) in units of MeV, i.e. the level density is 0.67 (RMF), 0.81 (RMF+PCV) 
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FIG. 6: The typical strongly fragmented states far from the Fermi surface in the odd mass nuclei 
surrounding 208 Pb calculated within RMF with allowance for the particle-vibration coupling. 

and 0.94 (EXP) in units of MeV -1 . This corresponds to an effective mass m* /m w 0.71 for 
the RMF and m*/m ~ 0.85 for the RMF+PCV calculations. We observe that the values for 
the effective masses of the protons are slightly smaller than those of the neutrons. 

28[ have discussed in great detail the concept of the effective mass 



Jaminon and Mahaux 
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in the case of RMF theory. On one side one has the well known Dirac mass 

m D = m + t s (r), (45) 

which is determined by the scalar field S s . Since we do not use an iso- vector scalar field for the 
present parameter set NL3 the Dirac mass is in these calculations identical for protons and 
neutrons. However, this quantitiy should not be compared with the effective mass determined 
empirically form a non-relativistic analysis of scattering data and of bound states. From a 
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FIG. 7: Neutron single-particle states in Pb 208 : the pure RMF spectrum (left column), the levels 
computed within RMF with allowance for the particle-vibration coupling (center) and the experi- 
mental spectrum (right). 

non-relativistic approximation of the Dirac equations one finds that the mass 

m e ff = m - S (46) 

should be used for this purpose. Here S is the time-like component of the Lorentz vector 
field determined by the exchange of to- and p-mesons. 

In symmetric nuclear matter we find for NL3: mo/m = 0.60 and m e ff/m = 0.67. The latter 
value is smaller then the values m*/m pa 0.71 for protons and m* jm pa 0.76 for neutrons 
deduced from the calculated spectrum around the Fermi surface in simple RMF theory in 
Figs. E|and|Bl Following similar arguments we would obtain for RMF+PVC calculations an 
average effective mass of 0.89. This is obviously still too low as compared to the experimental 
value. 

From the other hand, around the Fermi surface where relativistic kinematic effects are not 
significant our RMF+PVC spectrum can be characterized by the effective mass deduced 
from the Schrodinger equation which is a non-relativistic limit of the Dirac equation (jSJ. 
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FIG. 8: The same as in Fig. but for proton single-particle states. 

In this approximation one can calculate the state-dependent E-mass m/m RMF which is the 
inverted spectroscopic factor of the dominant level A: 

^=Kr- (47) 

For the calculated RMF+PVC spectrum the averaged E-masses are 1.26 for neutrons and 
1.41 for protons if one takes into account all the states given in the Table ITTT1 with spectro- 
scopic factors more then 0.5, i. e. good single-particle states. Thus, the energy dependence 
of the mass operator increases the RMF neutron and proton effective masses up to the values 
0.96 and 1.0, respectively. 

C. Comparison with other approaches 

Although the problem of particle-vibration coupling in nuclei has a long history and it was 
considered in a number of works, most of them are based on a non-relativistic treatment of the 
nuclear manv-body p r ob,„, 0„l y in a relatively recent investigation in Re, Q a correction 
of the RMF single-particle spectrum was undertaken in a phenomenological way assuming a 
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linear dependence of the mass operator near the Fermi surface. The corresponding coupling 
constants were determined by a fit to nuclear ground state properties. Despite the fact 
that the present approach is fully microscopic without any additional parameter adjusted to 
experiment, we find good agreement with the results of Ref. j3] for the spectrum of 208 Pb. 
The shift caused by the phenomenological particle-vibrational coupling in Ref. 2t| is only 
slightly larger than in the present investigation. 

Non-relativistic microscopic investigations of particle-vibrational coupling can be divided 



into two major groups. The first group |2i 



H uses a phenomenological single- 



particle input to reproduce the experimental spectrum and has therefore to exclude the 
contribution of the particle- vibration coupling from the full mass operator to find the 'bare' 
spectrum. Usually these older approaches take into consideration only a relatively small 
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number of collective low-lying phonons and use a particle-vibration coupling model 
This restriction to only low-lying modes produces shifts less then 1 MeV. However, as it was 
shown in Ref. [sj], enlarging of the phonon space with high- lying vibrations leads to very 
strong shifts of the single-particle levels up to 4 MeV, and no saturation is observed with 
respect to the dimension of the phonon space. 
The second group of approaches (see, for instance Refs. 

HQ) 

starts from a self-consistent 

Hartree-Fock description and applies perturbation theory to calculate the particle- vibration 

contribution to the full mass operator. In such self-consistent methods it is more justified 

I — I 

to enlarge the phonon space. It was shown, for instance, in |33| that the contribution of 
the isovector modes is noticeably smaller than the isoscalar ones. The detailed investigation 
of the relative importance of the high multipole states was performed in 3^. Because of 
the larger phonon space the typical shifts of the single-particle levels in 208 Pb are about 1-2 
MeV. 

As for the spectroscopic factors, all the approaches predict similar values because these 
factors are not very sensitive to the details of the calculational schemes. 
Thus one can see that the results of the present work are in a good agreement with the 
results of earlier approaches. 
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IV. SUMMARY 



The problem of the particle-vibration coupling is considered on the foundation of the rel- 
ativists mean field approach. The Dyson equation for the exact single-particle Green's 
function is solved in the Dirac basis by taking into account the energy-dependent part of 
the fully relativistic mass operator. This energy- dependent part is treated in terms of the 
particle-vibration coupling model that has been applied for the relativistic approach. 
The particle-phonon coupling amplitudes have been computed within self-consistent RRPA 
using the parameter set NL3 for the Lagrangian. A rather large number of collective vibra- 
tional modes has been taken into account. Relativistic effects of the Dirac sea on the mass 
operator have been analyzed in the usual no-sea approximation. They have been found to 
be small as compared to the contribution of the states with positive energy. Nontheless the 
the Dirac sea contributions are crucial for the description of the RRPA vibrations (Ref. [lfil ]) 
and are therefore fully taken into account. 

Noticeable increase of the single-particle level density near the Fermi surface relatively the 
pure RMF spectrum is obtained for 208 Pb that improves the agreement of the single-particle 
level scheme with experimental data considerably. For four odd mass nuclei surrounding 
208 Pb the distribution of the single-particle strength has been calculated and compared with 
experiment as well as with the results obtained within several non-relativistic approaches. 
The major result of the present work is a consistent description of nuclear many-body dy- 
namics including complex configurations within an approach which is (i) fully self-consistent, 
(ii) based on relativistic dynamics, (iii) universally valid for nuclei all over the periodic table, 
and (iv) based on a modern covariant density functional, which has been applied with great 
success of many nuclear properties all over the periodic table. Complex configurations play 
an important role in our understanding of the dynamics of the nuclear many-body problem. 
Here we have discussed the single particle motion and its coupling to collective vibrations. 
Such configurations are also of great importance in the description of damping phenomena 
in even-even nuclei. Thus, it is also interesting to investigate how the coupling to vibrational 
states produces the spreading width within an extended relativistic RPA approach. Work 
in this direction is in progress. 
ACKNOWLEDGEMENTS 

We are thankful to Dr. V. I. Tselyaev and D. Pena Arteaga for useful comments and 



25 



discussion. This work has been supported in part by the Bundesministerium fur Bildung und 
Forschung under project 06 MT 193. E. L. acknowledges the support from the Alexander von 
Humboldt-Stiftung and the assistance and hospitality provided by the Physics Department 
of TU-Munchen. 



[1] M. Bender, P.-H. Heenen, and P.-G. Reinhard, Rev. Mod. Phys. 75, 121 (2003). 
[2] P. Ring, Prog. Part. Nucl. Phys. 37, 193 (1996). 

[3] D. Vretenar, A. V. Afanasjev, G. A. Lalazissis, and P. Ring, Phys. Rep. 409, 101 (2005). 

[4] T. Gonzales-Llarena, J. L. Egido, G. A. Lalazissis, and P. Ring, Phys. Lett. B379, 13 (1996). 

[5] J. Boguta and A. R. Bodmer, Nucl. Phys. A292, 413 (1977). 

[6] Y. K. Gambhir, P. Ring, and A. Thimet, Ann. Phys. (N.Y.) 198, 132 (1990). 

[7] G. Lalazissis, D. Vretenar, and P. Ring, Eur. Phys. J. A22, 37 (2004). 

[8] G. A. Lalazissis, M. M. Sharma, P. Ring, and Y. K. Gambhir, Nucl. Phys. A608, 202 (1996). 

[9] J. Meng and P. Ring, Phys. Rev. Lett. 77, 3963 (1996). 
[10] G. A. Lalazissis, D. Vretenar, and P. Ring, Phys. Rev. C 69, 017301 (2004). 
[11] G. A. Lalazissis, D. Vretenar, and P. Ring, Nucl. Phys. A650, 133 (1999). 
[12] A. V. Afanasjev, J. Konig, P. Ring, L. M. Robledo and J. L. Egido, Phys. Rev. C 62, 054306 
(2000). 

[13] A. V. Afanasjev, P. Ring, and J. K6nig, Nucl. Phys. A676, 196 (2000). 

[14] Z.-Y. Ma, A. Wandelt, N. Van Giai, D. Vretenar, P. Ring, and L.-G. Cao, Nucl. Phys. A703, 
222 (2002). 

[15] D. Vretenar, H. Berghammer, and P. Ring, Nucl. Phys. A581, 679 (1995). 
[16] P. Ring, Z.-Y. Ma, N. Van Giai, D. Vretenar, A. Wandelt, and L.-G. Cao, Nucl. Phys. A694, 
249 (2001). 

[17] N. Paar, P. Ring, T. Niksic, and D. Vretenar, Phys. Rev. C 67, 034312 (2003). 
[18] A. Ansari, Phys. Lett. B623, 37 (2005). 

[19] P. Ring and P. Schuck, The nuclear many-body problem (Springer, Heidelberg, 1980). 
[20] J. P. Jeukenne, A. Lejeunne and C. Mahaux, Phys. Rep. 25C, 83 (1976). 
[21] A. B. Migdal, Theory of Finite Fermi Systems and Applications to Atomic Nuclei (Interscience, 
New York, 1967). 



26 



[22] A. Bohr and B. Mottelson, Nuclear Structure, Vol. 2 (Benjamin, New York, 1975). 

[23] S. P. Kamerdzhiev and V. I. Tselyaev, Yad. Fiz. 44, 336 (1986) [Sov. J. Nucl. Phys. 44, 214 
(1986)]. 

[24] P. Ring and E. Werner, Nucl. Phys. A211, 198 (1973). 

[25] G. A. Lalazissis, J. K6nig, and P. Ring, Phys. Rev. C 55, 540 (1997). 

[26] A. Avdeenkov and S. Kamerdzhiev, Phys. Lett. B499, 423 (1999). 

[27] M. Jaminon and C. Mahaux, Phys. Rev. C 40, 354 (1989). 

[28] M. Jaminon and C. Mahaux, Phys. Rev. C 41, 697 (1990). 

[29] D. Vretenar, T. Niksic, and P. Ring, Phys. Rev. C 65, 024321 (2002). 

[30] A. P. Platonov, Yad. Fiz. 34, 612 (1981) [Sov. J. Nucl. Phys. 34, 342 (1981)]. 

[31] I. Hamamoto and P. J. Siemens, Nucl. Phys. A269, 199 (1976). 

[32] R. P. J. Perazzo, S. L. Reich, and H. M. Sofia, Nucl. Phys. A339, 23 (1980). 

[33] V. Bernard and N. Van Giai, Nucl. Phys. A348, 75 (1980). 



27 



